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Abstract 

A new approach to multi-dimensional quantum scattering by the infinite 
order discrete variable representation is presented. Determining the expansion 
coefficients of the wave function at the asymptotic regions by the solution 
of the differential Schrodinger equation, we reduce an infinite set of linear 
equations to a finite one. Application to the benchmark collinear H + H2 —* 
H2 + H reaction is shown to yield precise reaction probabilities. 
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One of the most common approach to the solution of quantum scattering problem is 
the application of square-integrable (L 2 ) basis functions fTj-Q. Of the L 2 basis methods the 
discrete variable representation(DVR) method || is proven to be highly successful 0-0. 
The DVR is a grid-point representation in which the potential energy is diagonal and the 
kinetic energy is a sum of one- dimensional matrices. Hence the Hamiltonian matrix is 
extremely sparse, which means that iterative linear algebra methods [HJ can deal efficiently 
with extremely large systems. 

Recently, Eisenberg et. al. [11] has developed an infinite order DVR method for one- 



dimensional quantum reactive scattering problems. They expanded the wave function in 
terms of an infinite set of L 2 basis functions satisfying the conditions of DVR P-^|T2"| . The 



matrix related to the resulting set of algebraic equations, though infinite, has been shown to 



have the structure of a Toeplitz matrix . Using the analytical properties of the Toeplitz 
matrix, they reduced the infinite set of algebraic equations to a finite one and obtained very 
accurate results for one-dimensional potential scattering. However, their method is not likely 
to be extended to multi-dimensional systems due to the failure of the Toeplitz properties. 

In this Letter, we show that the infinite order DVR can be properly converted into a 
finite DVR even for multi-dimensional reactive systems. Instead of investigating the ana- 
lytical properties of the Toeplitz matrix, we use the solutions of the differential Schrodinger 
equations at the asymptotic regions. 

The full scattering wave function \l/ n is decomposed of the incoming distorted wave 0~ 
and the outgoing wave xt, j 

= -0- + (1) 

where n is a channel index, i.e., a superindex over the arrangement and rovibrational in- 
dices. The distorted wave 6~ is any regular scattering solution corresponding to a simple 
nonreactive (i.e., arrangement conserving) Hamiltonian, H 0n , in the asymptotic region in 
the arrangement channel n: 

(H 0n - E)6- = 0, (2) 



while \l/ n is the solution of 

(H-E)* n = 0, (3) 

where H is the full Hamiltonian. For the sake of later convenience, we impose the "totally 
incoming" boundary condition on 9~ as 

0-(r,i2)~— — u^r) (4) 

sJVn 

for large R. Here r denotes all possible internal coordinates of the system, thus excluding the 
channel radius (relative translation) R; {u n (r)}, the channel eigenfunctions; v n the relative 
velocity for channel n, and k n the corresponding wave vector. The equation for xt i s then, 
from (l)-(3), 

{E-H)xt = {H-E)9-. (5) 

Xn must obey "totally outgoing" boundary condition, i.e., the only incoming wave part in 
the full wave function is due to 9~. 

Using the infinite order uniform DVR ||, (5) is transformed into an infinite set of cou- 
pled linear algebraic equations. To do this, we introduce the following convenient sets of 
DVR basis functions Q, 7Z , Vq, R and V as shown in Fig. 1: Q be the set of N DVR 
basis functions represented by the N grid points in the reactive (strong) interaction region 
correspond; TZq = {H^l), ...,lZo(N r )}, the set of N r DVR basis functions in the reactant 
asymptotic (relatively weak interaction) region nearest to the region of reactive interaction; 
Vo = {Pq(1), ...,Vo{N p )}, the set of N p DVR basis functions in the product asymptotic re- 
gions nearest to the region of reactive interaction; TZ, the infinite set of DVR basis functions 
in the reactant asymptotic region except the functions in TZq, and V the infinite set of DVR 
basis functions in the product asymptotic regions except the functions in Vq. Here N r (N p ) 
is the number of open channels in the reactant (product) arrangement. Using these sets, we 
can rewrite (5) as a set of coupled algebraic equations, 

A non {TZ\xt) + ^WoMxt) + An oQ (Q\xt) + ^n Q v (V \xt) + ^W<P|x£> = (^o|# " E\®nX ( 6 ) 
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A Q n{K\xt) + A Q K {Ko\xt) + Ma{Q\xt) + A Qv (^o\xt) + A Qv(V\xt) = (Q\H - E\9^), (7) 
Av n(K\xt) + A Vo iio(K \xt) + A r Q(Q\xt) + A v v (V \xt) + A PoP (V\xt) = W H ~ E \ d n^ ( 8 ) 

where 

A ij = {E-V j )5 ij -T ij , (9) 

Tij is the kinetic energy matrix element ||, which is analytically obtained, connecting the 
DVR grid points % E TZ + Q + V and j E TZ + 7Z + Q + V + V, Vj is the potential 
energy at the DVR grid point j, and we omit the summation over the index j such that, 

e.g., A non (K\xZ r ) means E j& n+n Q +Q+r +v A ijU\xt) for i e n o- 

In the above we did not write explicitly the similar equations corresponding to the TZ 
and V component since we do not use them in order to eliminate (7Z\Xn r ) an d (P\xt, r ) by 
searching for an analytical property as in the case of Toeplitz. Instead, for the asymptotic 
regions, we introduce, 

N p 

CP\xt) = ~ E (^lt>V. (io) 

n p =l 

N r 

(K\xt) = ~ E (KK r )S K nr, (11) 

n' r =l 

where 9+ is the regular "totally outgoing" wave which satisfies (2), i.e., (V\9+) = (V\9~)*, 
and S is the scattering matrix to be determined. (10) and (11) can be rewritten as 

N p N p 

<P\xt) = E EWi> (t v Y . (v ®\xt)> (i 2 ) 

. * \ / 'lip I 

n p =l 1=1 

N r N r , 

(K\xt) = E EWJ.) (T-*)- (K (i)\xt), (13) 



n'=l i=l 



where 



TL = m)K) (1 4 ) 



and 



T^ = (7Z (z)\9t). (15) 



The above equations (12)-(15) are the key in this Letter. Substituting (V\Xn r ) an d (T^\Xn r ) 
in (12) and (13) into (6) through (8), we obtain the following set of N + N r + N p linear 
equations: 

(n \H-E\e- r ), (16) 

(Q\H-E\9- r ), (17) 
(V \H-E\9- r ), (18) 

(19) 

(20) 

where % is a DVR grid point in the sets 7£o, <2, and V§. Note that e^- is nonvanishing only 
for j in the sets 71$ and Vq. Applying the complex conjugate of (12)-(15), (19) and (20), 
the right hand sides of expressions (16) through (18) can be expressed as 

(i\H - E\e- r ) = - (a + e*) iUo (n \e- r ). (21) 

Hence, (16)-(18) can be expressed in the form of an (N + N r + N p ) x (N + N r + N p ) matrix 
equation 

A- X + =-A*.0~, (22) 

where A = A + e. Let G + be the left-inverse of A, i.e., G + ■ A = 1. Then it can be shown 
that 

Xt = ~K + G + • (e - e*) ■ 0~ r . (23) 

Finally we obtain the full scattering wave function and the scattering matrix, 

vtV = G + .(6-e*).0-, (24) 
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(e + A) no KMxt) + A n oQ (Q\xt) + (e + A )n v (r \xt) 
(e + A) Qno (n \xt) + A Q Q(Q\xt) + ( e + A )aPo(Po\xt) 
(e + A) V()7lo (n \xt} + a Voq(Q\xZ) + ( e + ^Wo(^o|xt> 

where e is the matrix of which the non-vanishing elements are 



jen n' r 



-i 
no 



and 

^ = EEA,0K)(3 VP ) B "! 



ripVo 



and 

Sn pnr = E [G + • (6* - 6)] (A" ) (25) 

For the calculation of e, we need the overlaps CJZ\0^, ), (TZo^, ), (^1^)) an< ^ (^ol^n )■ Using 
the property of the DVR basis functions, 

{i\%) = e+{r il R i )Vfi 1 (26) 

where r { and Ri are respectively the values of r and at the DVR grid point i, and w is the 
quadrature weight |5|-^,12|, the above overlaps can be easily obtained by numerical integra- 



tion of the differential Schodinger equation (2) corresponding to the nonreactive scattering 
region discussed earlier. 

Thus far we have shown that the derived scattering matrix (25) is completely general 
enough to deal with the fully three dimensional scattering processes. For the sake of compar- 
ison with other accurate calculations, we apply the method to the collinear H + H2 — > H2 + 
H reaction in the energy range of 0.4—1.0 eV. The coordinate axes for the DVR are chosen 



to be the normal ones of the transition state [16[ as shown in Fig. 1. We use the LSTH 



potential energy surface [14|. The distorted wave functions, 9+, are calculated by integrating 
the coupled differential equations (2) inward from R asym p(> 20.0 a.u.) to R max (> 6.0 a.u.) 
for each arrangement, and then are interpolated at the DVR grid points in the asymptotic 
regions to obtain the overlaps (Tl\9^, ), (TZq\9^, ), etc. Ei—0. 78629 eV is the threshold of the 
vibrational quantum states with n = 1. Hence there are two linearly independent distorted 
wave functions 9q and 9f at the energies above E\. 

Fig. 2 illustrates the convergence of Pq(E) and P\(E) at E = 0.9678 eV, where P nr (E) 
is the total reaction probability from a given reactant molecular state n r to any product 
molecular states at energy E, i.e., 

P nr (E)= £ I^J 2 , (27) 

n p =l 

as the size of the region of reactive interaction, R ma x, is increased, for different values of 
the grid constant ub [pfL6|. The percent error is calculated by comparison to the results of 



Colbert and Miller ||. The convergence for E = 0.9678 with respect to Ub for Rmax = 8.0 
is also shown in Fig. 3. The total reaction probabilities, with R max = 8.0 and = 3.4, 
are converged to within 0.05% at this energy,. Fig. 4 shows excellent agreement of the 



converged results for Pq{E) and Pi(E) with the exact results of Bondi and Connor [HJ over 
a range of total energies including the resonance energy near 0.9 eV. 

In summary, the infinite order DVR theory for multi-dimensional reactive scattering has 
been presented. By applying it to the collinear H + H2 — > H2 + H reaction, we obtained 
precise reaction probabilities comparable with the other results obtained from the direct 
integration of the differential Schrodinger equation (3) ||15|| . The advantages of the present 
theory are as follows: 1. iterative linear algebra methods are readily applicable, 2. no 
numerical integrals are needed for calculating the matrix elements, 3. as can be seen from 
(19), (20), (24) and (25), no regularization techniques |3||| for the distorted wave functions 
are necessary since no informations on those functions in the reactive (strong) interaction 
region are required although they are set to be regular at R = in the beginning (1), 4. the 
effective Green's operator G + in the reactive interaction region is obtained as a by-product, 
5. the formalism is completely general and can be applied to any non-collinear scattering 
problems including collinear systems. 

We would like to thank D.G. Choi for helpful discussions and W.H. Thompson for pro- 
viding information on the calculational results in references [15] and ]TjJ. This work was 
partially supported by the Center for Molecular Science at KAIST, and the numerical cal- 
culations were performed on the computer CRAY-C90 at SERF 
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FIGURES 

FIG. 1. Schematic diagram showing the sets 1Z (plus), TZq (diamond), Q (dot), Vq (square), 
and V (cross) of the DVR basis functions, which are represented by the DVR grid points, for the 
collinear H + H2 — ► H2 + H reactive scattering problem. Three straight lines are drawn to divide 
the sets 1Z, Q, and V from each others. The potential energy contours and the (x, y) — coordinate 
system for the DVR are also depicted. 

FIG. 2. Relative percent error of the total reaction probabilities Pq(E) (solid lines) and P\{E) 
(dashed lines) at E = 0.9678 eV as a function of R max for different values of ub\ ng = 3.4 (thick 
lines); and ub = 1.7 (thin lines). 

FIG. 3. Relative percent error of the total reaction probabilities Pq{E) (solid lines) and P\{E) 
(dashed lines) at E = 0.9678 eV as a function of ng for R max = 8.0. 

FIG. 4. Total reaction probabilities (a) Pq(E) and (b) Pi(E) as a function of energy computed 
using Rmax = 8.0 a.u. and ng = 3.4 (circle) compared to the results of Bondi and Connor (solid 
line) (ref. [15]) and Thompson and Miller (triangle) (ref. [16]). 
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Fig. 1 




Fig. 2 
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Fig. 3 
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Fig. 4(b) 
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